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Abstract 

o 

OA We provide a complete characterization of the class of one-dimensional time- homogeneous 
diffusions consistent with a given law at an exponentially distributed time using classical 
results in diffusion theory. To illustrate we characterize the class of diffusions with the same 
distribution as Brownian motion at an exponentially distributed time. 

i— - 1 1 Introduction 

The aim of this article is to characterize the class of one-dimensional time-homogeneous diffu- 



sions with a given law at an exponentially distributed time. We show, for instance, that there 

B 



> 



is a one-parameter family of diffusion processes started at with the same law as Brownian 
motion at an exponentially distributed time. In general, given a probability distribution we 
find that consistent diffusions are parameterized by a choice of starting point and secondly by 
a choice of Wronskian. 



We use classical results due to Dynkin [4J and Salminen [IT] involving the /i-transform (or 
Doob's /i-transform) of a diffusion to provide necessary and sufficient conditions for a diffusion 
to have a given distribution at a random time. Previously, Cox, Hobson and Obloj [3] proved the 
existence of consistent diffusions when the first moment is finite. We recover the construction 
in [3j as a canonical choice from the class of consistent diffusions. 

The problem of constructing diffusions with a given law at a random time can be motivated 
as an inverse problem in finance. Suppose we are given European call option prices with a fixed 
expiry time for a continuum of strikes. It is a natural inverse problem to ask whether there exist 
models for the asset price process consistent with these call prices and desirable properties. This 
article is related to a particular solution of the inverse problem proposed by Carr and Cousot 
[2]. The idea in [2] is to construct gamma-subordinated martingale diffusions to fit call prices 
after scaling the gamma-subordinator to be exponentially distributed at the expiry time. 

It is well known (see Breeden and Litzenberger [Tj ) that knowledge of call prices is equivalent 
to knowledge of the marginal law of the underlying asset. The problem of constructing a gamma- 
subordinated martingale diffusion consistent with call options at a single maturity is therefore 
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a version of the more general problem considered here. In this article we do not assume that 
the consistent process is a martingale, nor that the target distribution is centred and we recover 
a parameterized class of consistent diffusions rather than a particular consistent martingale 
diffusion. 

The analogue problem of constructing diffusions with a given distribution at a deterministic 
time is considered by Ekstrom, et. al. in [5i|. For a more general view on inverse problems see 
also Hobson [7J. 

Finally, this article is related to the inverse problem of constructing diffusions consistent 
with prices for perpetual American options or, more generally, with given value functions for 
perpetual horizon stopping problems, see Hobson and Klimmek [6j. As in this article, the 
underlying key idea in [6] is to construct consistent diffusions through the speed measure via 
the eigenfunctions of the diffusion. 



2 Generalized diffusions and the /i-transform 

Let I C M be a finite or infinite interval with a left endpoint a and right endpoint b. Let m be 
a non-negative, non-zero Borel measure on R with / = supp(m). Let s : / — > M. be a strictly 
increasing and continuous function. Let B = (Bt)t>o be a Brownian motion started at -Bo = 
s(Xq) supported on a a filtration ¥ B = (J 7 ^) u >o with local time process {L z u ]u > 0,z € M.}. 
Define T to be the continuous, increasing, additive functional 

T u = Llm(dz), 

JR 

and define its right-continuous inverse by 

A t = inf{u : T u > t}. 

If Xt = s^ 1 (B(At)) then X = (Xt)t>o is a one-dimensional regular diffusion with speed measure 
m and scale function s and Xt £ I almost surely for alH > 0. 
Let H x = inf{n : X u = x}. Then for A > (see e.g. [IT]). 



My) x ~ y > 



Ux,y) = E x [e- XH y] = \ £ A fj x z y (2.1) 



where ip\ and <j)\ are respectively a strictly increasing and a strictly decreasing solution to the 
differential equation 

1 d d „ , „ , 
r/ = V- ( 2 - 2 ) 

The two solutions are linearly independent with Wronskian W\ = <p'x4>\ ~ ftx^x > 0- Recall 
that if a diffusion X = (Xt)t>o is in natural scale, then the Wronskian W\ is a constant. In 



the smooth case, when m has a density v so that m{dx) = v{x)dx and s" is continuous, (2.2) is 
equivalent to 

l -e 2 {x)f\x) + a(x)f'(x) = Xf(x), (2.3) 



where 



v (x) = a~ 2 {x)e Mi > x \ s'{x) = e- M{x \ M{x) = [" 2a~ 2 (z)a(z)dz. 

Jo- 
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We will call the solutions to ( |2.2[ ) the A-eigenfunctions of the diffusion. We will scale the 
A-eigenfunctions so that <p\(Xo) = <fi\(Xo) = 1. 

The A-eigenfunctions are well known to be A-excessive. We recall that a Borel- measurable 
function h : I — > R + is A-excessive if for all x £ I and t > 0, Kx[ e ~ Xt h(Xt)] < h{x) and if 
E x [e~ Xt h(X t )} -4 pointwise as i ->■ 0. 

Definition 2.1. Let /i be a A-excessive function. The h-transform of a diffusion X = (Xt)t>o 
is the diffusion X h = (X[ ;)t>o with transition function 

P h (t;x,A) = -L /" e - A *p(t ; x,y)/i(y)m(dy), 



where p is the transition density of X. 

By the following result due to Dynkin [3] (see also Salminen (3.1)), any diffusion X can 
be transformed into a diffusion with a given marginal at an exponential killing time. Fix A > 
and let X = (Xt)t>o be a diffusion with A-eigenfunctions ip\ and <p\. Let T be an exponentially 
distributed random variable with parameter A, independent of X. 

Theorem 2.2. Given a probability measure \jl on [a,b] let 

h(x) = [ jT^Kdy)- (2.4) 
J[a,b] €\{Xo,y) 

Then P(X^ £ dx) = fi(dx). Conversely, let h be a X-excessive function with h(Xo) = 1 and let 



Jx{dy) = F{X^, G dx). Then h has the representation (2.4) with fx = 7^. 

The measure 7^ in (2.4) is called the representing measure for h. It follows from Theorem 



2.2 that we can start with any diffusion X on [a, b] and construct a killed diffusion with a given 
representing measure via an /i-transform. Thus, in principle, since the representing measure co- 
incides with the law of Xj,, Dynkin's result solves the inverse problem of constructing diffusions 
with a given marginal at an exponentially distributed (killing) time. 

We will build on this observation to recover consistent diffusions using a characterization of 
a representing measure in terms of the A-eigenfunctions. 

3 Characterizing consistent diffusions 

Without loss of generality, we will restrict the inverse problem to the class of diffusions in natural 
scale. If it is possible to construct a diffusion in natural scale consistent with a given law on 
[a, b] at an exponential time then it follows that we can also construct consistent diffusions 
with non-trivial scale functions: Given a marginal density \i on M, define a measure v via 
v(B) = fj,(s~ 1 (B)) for Borel sets B C M, where s : [a,b] — > M is an arbitrary scale function. 
Then if a diffusion Y = (Yt)t>o m natural scale is consistent with v, the diffusion X = s(Y) is 
consistent with [i. 

Recall the definition of the /i-transform and observe that h = 1 is a A-excessive function 
for any A > 0. We will call the /i-transform corresponding to h = 1 the A-transform. The 
A-transform of X is equivalent to X up to the exponential time T ~ Exp(X) at which time X 1 
is killed, while X remains on the state space I. Thus 

n- 
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where A is the grave state of the killed diffusion X . Note that the transition density of X 1 
is given by q(t;x,y) = e~ xt p(t; x,y), where p(t;x,y) is the transition density of X. Other 
fundamental quantities are related similarly, for instance Ex [e~ XHx ] = Fx {H x < T) = 
PxoiX 1 reaches x). 

We restate our inverse problem as follows. Given a probability measure [i on [a, b], construct 
a diffusion X = (Xt)t>o such that for all x G [a, b] 

&(X, *W (3-D 



0,6] t\( x o,y)' 



whence by Theorem 2.2 X T ~ /j,. Since Xt = X T ~ //, the idea is to construct the class of 
consistent diffusions via the A-eigenfunctions for which (3.1) holds. 

The following result is an elementary case (h = 1) of Proposition (3.3) in Salminen 

Proposition 3.1. Given a diffusion X, the representing measure 7 = 7^ is given by 

<d (x-) 

j([a,x)) = ^-—>- fora<x<X , (3.2) 

J((x,b}) = ~^ X+) forX <x<b, (3.3) 



where (fx (<p\) are the increasing (decreasing) X-eigenf unctions of X and W\ is the Wronskian. 
Remark 3.2. Suppose that a is accessible and that Xq = a. Then the representing measure 



for h = 1 is given by j((x,b\) = ^^} X+ ^ for a <x <b and similarly if Xq = b, where b is 



accessible. 



The characterization of the representing measure given by Proposition 3.1 will be used 
to arrive at our main result. Suppose we are given a probability measure fi on [a,b\. Let 

U ^ x ) = I[a,b] \ x - yH*/), CVO) = f[afi](v - X ) + Kdy) and P^x) = j [aM (x - y)+{i{dy). Let 
X = (Xt)t>o be a one-dimensional diffusion in natural scale and let T be an independent 
exponentially distributed random variable with parameter A > 0. 

Theorem 3.3. Suppose Xq G (a, b). Then Xt ~ /i if and only if the speed measure of X 
satisfies 



m(dx) 



1 /i(dx) <^ ~Y~ 

2AP^(x)-P M (X )+l/H/ A < X - 

J M^) v <r ^ ^ h 

2A C Al (x)-C M (Xo)+l/W A ^O^X <0. 



where W\ > is the Wronskian of X. 



Proof. Suppose first that Xt ~ fJ.- Then since Xt = X T , by Theorem 2.2 \x is the representing 
measure for h = 1. Differentiating both sides of (3.2) we find that for all points x such that 
a < x < Xq and which are not atoms of /i, 

Kdx) = ^(fx(x)dx. 

(If ji has an atom at x then /x({x}) = yy-(ip'^(x+) — ip'^(x—)). The case x > Xq is similar, with 
4>\ replacing <p\. 
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On the other hand, integrating the two sides of (3.2) we have 



P^{x) + k l 



W x 

4>\{x) 



for x < Xq, 
for x > Xq, 



for constants ki,k 2 G M. Now using the fact that 99^(^0) = ^a(^o) = 1 we find that k\ = 
1/W\ — P(Xq) and k 2 = 1/W\ — C(Xq). Since tp\ and <t>x are the A-eigenfunctions for X and 



solutions to (2.2), the speed measure of X satisfies 

m{dx) = 



1 <p"{x)dx 
2A Vx (x) 
1 <j>'{{x)dx 



2A </> A (:r) 

Substituting for ipx and </>a we thus have 

C 1 fi{dx) 



m{dx) 



2A P M (x)+fci 

1 



a < x < Xq 
Xq < x < b. 

a < x < Xo 
X < x < b 



I 2AC fl (z)+fc 2 

as required. 

Conversely suppose that X has the given speed measure on (a,b). Define a function rj : 
[a,b] — y M + as follows. Let W\ > be the Wronskian associated with X and set 



rj{x) 



Wx(P^x) - P^(xq)) + 1 a < x < Xq 
Wx(C^(x) - C^Xq)) + 1 Xq < X < b. 



Then n solves (2.2) on the domain (a, b) and we therefore have 

V(x) = 



<P\(x) a < x < Xq 
4>x{x) Xq < x < b. 



By Proposition 3.1 the representing measure for h = 1 is given by 



7([o,a?)) 
7 ((x,6]) 
and it follows that Xt ~ fi- 



rj'(x—) 

W x 
—rf(x+) 



fi([a,x)) for a < x < Xq 
■ fi((x, b]) for Xq < x < b, 



□ 



Remark 3.4. If X is started at an accessible end-point, a say, then Xt ~ H if and only if for 

all x G [a, b), m(dx) = jx c (x)-'c d (a)+i/w x ' ^ e case ^° = ^ wnere & is accessible is analogous. 
Compare Remark |3.2| 



We have the following interpretation for the Wronskian. 
Corollary 3.5. If Xt ~ n then the Wronskian satisfies 

Wx m(dz) 



2A 



fi(dz) 



z=X 



Intuition for Corollary 



3.5 



is provided by the fact that 2/Wx = Exot^^y] ( see Lemma VI. 
54.1 in Rogers and Williams jlOj). 
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4 The Wronskian and the martingale property 



By Theorem 3.3 the class of diffusions with a given starting point and a given law at an 
exponentially distributed time is parameterized by a choice of Wronskian. We will now see that 
when the first moment of the target law is finite there exists a unique consistent diffusion with 
the property that it is a martingale away from the boundary. 

We will suppose for the remainder of this section that Jj afe j \x\fi(dx) < oo. Let = 
fy ab ^x[i(dx) and let us fix the starting point, Xq = x^. Let r = inf{i > : X t $ int(I)}. 
Then X T = (Xt/\ T )t>o is a local martingale. In the following we will say that X is a martingale 
diffusion whenever X T is a martingale. 



It follows from Theorem 3.3 that X is a diffusion with law fi at an exponentially distributed 
i only 

m(dx) 



time T ~ Exp{\) if and only if 

I [i(dx 



A U, L (x) -\x- X \ - 2C /1 (X„) + 2/W x 
for a < x < b. Given this formula for the consistent speed measures, the most natural choice of 



W\ is perhaps W\ = 1/C fl (x fl ). Indeed, will see in Theorem 4.1 below that this choice recovers 



the unique consistent martingale diffusion (it also recovers the construction in [3]). 

Theorem 4.1. Suppose Xq = x^ and a = — oo or b = oo. Then X is a martingale diffusion 
consistent with /i if and only ifW\ = 1/C^(x^). 

The author would like to thank David Hobson for providing the proof that f°° dx = oo 

which is used below. 

Proof. We suppose b = oo, the case a = oo is analogous. Since m is positive, W\ > l/C^x^). 
Suppose W\ > 1/C^Xp) then 

1 n(dx) 

m[dx) - 



X U^(x) -\x - x^\+ c 



for some c > and lim — 7 - ) = 1/Ac. Thus f°° \x\m(dx) oc f°° \x\a(dx) < oo. It follows from 
ztoo n(dx) 

Theorem 1 in Kotani [9] that X is not a martingale diffusion. 

Conversely suppose that W\ = 1/C^(x^). We will show that J°° dx = oo. Write 



h(x) = For fixed y and x > y, let D(x) = E x exp (- f" y ^^dsj . Note that D(y) = 1 

and D is positive and decreasing. Let M t = exp ^— J * h ^ s ^ D(B t ). Then M = (Mt/\H y )t>o 
is a bounded martingale. In particular, by Ito's formula, ^D"(B S ) = h ^ s "> D(B S ), so that 
D(x) = ^H. It follows that lim D(x) = and that 



'* h(B s ) 



[Hy 

lim / - ds = oo 

xtoo Jo B s 
Bq=x 



6 



almost surely. Then we must have 



lim E 

xfoo 
Bo=x 

lim 

00 h(z) 

z 



H 



y h(B s ) 



ds 



x h(z) 



(z - y)dz + 



Hz) 



(x — y)dz 



(z - y)dz, 



and thus J°° h(z)dz = 00. By Theorem 1 in [9], X is a martingale diffusion. 



□ 



Remark 4.2. An alternative proof of Theorem |4.1| is available using a result in Hulley and 
Platen [8]. By Theorem 1.2 and Proposition 2.2 in [S], X is a martingale diffusion if and only 
if lim <j)\{x) = 0. Now recall that since X is consistent with /i we have 4>\(x) = W fi C IJi (x) — 

xfoo 

WaCJXo) + 1 for x > X . Clearly lim <f> x (x) = if and only if W„ = 1/C M (X ). 

xfoo 



5 Examples 

Example 5.1. Let B = (Bt)t>o be Brownian motion. Then we find that B? ~ fJ-x, where for 
x > 



(J>\({x, 00)) = /x A ((-oo, -x)) 



1 



2Ax 



Let us recover the class of consistent diffusions started at Xq = with the same law at an 
exponential time as Brownian motion. The consistent diffusions have speed measures mw{x) = 
vw(x)dx, where 

e -V2X\x\ 

u w (x) 



e -V2\\x\ - y f\/2 + 2\/W 

The choice W = 1/C(0) = 2\^2X corresponds to Brownian motion. Any choice of W G 
(0,1/C(0)) corresponds to a strict local martingale diffusion with the same marginal law. 




0.3 v w 



v w( x ) 



Figure 1: Plot of vw{x) for A = 1/2 and W G (0,2]. Note that v%{x) = 1 corresponds to 
Brownian motion which has Wronskian W = 2\/2A = 2. 



Example 5.2. Suppose that a = — 1, b = 1 and we wish to recover diffusions started at Xq = 
that are uniformly distributed at an exponential time. We find that the consistent diffusions are 
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parameterized by W S (0,4] with corresponding speed measures mw{dx) = vyy(x)dx given by 



\{x 2 + 2x + A/W) -l<x<0 
X(x 2 - 2x + A/W) 0<x<l 
oo otherwise. 



The canonical choice for W is 1/W = C(0) = 1/4. Since v±{— 1) = z^(l) = 00 the boundary 
points are inaccessible whence Xj. ~ U(— 1,1). For W G (0,4), the speed measure is finite on 
[— 1, 1] and the consistent diffusions reflect at the boundaries and Xjf ~ U[—l, 1]. 

Now suppose Xq = 1/2. Then 



\jv w (x) 



X(x 2 + 2x + 1/4 + 4/W) -l<x<l/2 
X(x 2 - 2x + 9/4 + A/W) l/2<x<l 
00 otherwise. 



2 
* 



Figure 2: Plot of fw(x) for A = 1/2 and W = 1 when = (solid line) and Xo = 1/2 (dashed 

w 



line), and Xq = —1 (alternating line). Note that v\y{Xo) = K = 1/2, see Corollary 



3.5 
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